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ABSTRACT 

We have used the Palomar 200" Adaptive Optics (AO) system to directly detect the astrometric 
brown dwarf GJ 802B reported by Pravdo et al. 2005. This observation is achieved with a novel 
combination of aperture masking interferometry and AO. The dynamical masses are 0.175±0.021 
M Q and 0.064±0.032 M Q for the primary and secondary respectively. The inferred absolute H band 
magnitude of GJ 802B is M#=12.8 resulting in a model-dependent T e fi of 1850 ± 50K and mass 
range of 0.057-0.074 M . 
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1. INTRODUCTION 

Binary stars provide a unique laboratory for the study 
of the physical properties of individual objects, and im- 
portant constraints on star formation and evolution. Al- 
though there are now a large number of objects known 
beyond the substellar limit, there are few in systems 
amenable to extracti on of dynamical meas urements of 
physical parameters l)Burgasser et al.l 120061) . This is a 
particularly acute issue in the case of substellar objects 
due to the degeneracy between age, mass and luminosity. 

We have selected targ ets from the STEPS survey 
ijPravdo fc Shaklanl Il996j) that show astrometric evi- 
dence of a low mass companion for imaging followup with 
Adaptive Optics (AO). In order to improve the sensitiv- 
ity to companions we have implemented a novel aper- 
ture masking interferometry technique in concert with 
the atmospheric turbulence correction afforded by AO. 
This technique has succeeded in directly detecting the 
low mass companion to the M5.5 dwarf GJ 802. 



2. OBSERVATIONS 



2.1. 



Aperture Masking Interferometry 

While the large gains of AO for high contras t imag - 
ing are widely reco g nised (e.g [Majdntos h et alJ |2003); 
ICarson et all pOOSTk iClose et all <|2005[) L nearly all AO 
searches for substellar companions have focussed on 
achieving very high dynamic range (> 10 4 ) at mod- 
erate separations. Several issues drive AO imaging to 
this parameter space. The interaction of the actuator 
count of AO systems and the limitations of coronag- 
raphy l|Sivaram akrishn an et alJl2001t iLlovd et alJl2001fl 
lead to both practical and fundamental limitations for 
high contrast imaging. Present and currently planned 
((Macintosh et alJ [2004) AO systems are focussed on 
achieving high contrast at radii of more than 4 X/D. 
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At closer separations, diffracted light is difficult to 
suppress with a coronagraph, and most importantly 
there is a large noise floor due to the presence of 
fluctuating speckles in the image iRacine et~al 119991: 
IFit.zgera.ld fc Orah ami 1207)3 ISoummer et tilEim. Fi- 
nally, it has proven to be remarkably difficult in practice 
to precisely calibrate the AO PSF. A variety of differ- 
ential imaging approches have been propos ed to circum- 
vent t he pr o blem of AO PSF calibration llMar ois e1j all 
1 20001 120041: iBrandner et all 120051: iMarois et all l2006t 
Sparks & Ford 120021) . These approaches rely on exploit- 
ing a differential signal in wavelength, polarization or 
sky rotation to improve the source extraction, but do 
not fundamentally address the issue of calibration of AO 
data. 



Telescope t AO System FHARQ 




Fig. 1. — Functional diagram of the aperture masking experi- 
ment. Although in practice the optical system is complex and in- 
volves many reflective optical elements, only the essential imaging 
properties are represented with lenses. After the telescope and AO 
system, the telescope pupil is reimaged at the PHARO Lyot stop. 
Under most common circumstances, the Lyot stop serves only as 
a cold baffle. In a coronagraphy mode, an undersized stop is used 
to block the light diffracted by the pupil edges. For the aperture 
masking interferometry mode, a pupil mask is placed in the Lyot 
stop to form an interferogram recorded at the focal plane. The im- 
ages shown are the full pupil AO images and 9-hole interferograms 
obtained for GJ 802 in 2004 September. 
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In light of these considerations, we have undertaken 
a novel experiment to achieve precision calibration of 
AO data, by marrying the sensitivity of AO observa- 
tions with the precision calibration afforded by interfer- 
ometry. The heritage of Non-Redundant Masking in- 
terferometry iTuthillet alJ E)0(1 iReadhead et alUloM 
iNakaiima et alJll989|) can be combi ned with the wave - 
front stabilization of adaptive optics ijTuthill et all2006|) . 
In practice, the optical implementation of this capabil- 
ity is relatively simple (see Figure QJ. We have used 
the Palomar 2 00" telescope with the PALAO adaptive 
optics system llTrov et all 120001) an d the PHARO in- 
frared camera llHavward et alJ 12001'!. PHARO was de- 
signed with coronagraphic capability in mind, so in- 
corporates a ten position Lyot wheel in the collimated 
beam at the internal position of the re-imaged telescope 
pupil. This wheel holds a variety of pupil stops to en- 
able the interchange of Lyot sto ps with various undersiz- 
ings QQppenheimer et"a l. 2000). We have installed non- 
rcdundant masks in the spare openings of the PHARO 
Lyot wheel. The 9-hole mask used in this work is opti- 
mized for broad band, faint targets with 50 cm diameter 
subapertures and 4.15 m longest baseline. The 9-hole 
mask transmits approximately 15% of the total light in- 
cident on the telescope pupil. 
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Fig. 2. — Extraction of closure phases from optical interfero- 
grams. Each "splodge" in the Fourier plane corresponds to the 
fringe formed by a pair of holes in the pupil mask. The splodges 
are point-symmetric since the interferogram is real (or equivalently 
each splodge appears for both the baseline formed by a pair of holes 
in both order A-B and B-A). Since the splodges map uniquely to a 
baseline, closure phase relations can be constructed that reject any 
residual phase errors, and therefore speckle noise. One such closure 
phase triangle (out of 84 possible for a 9 hole mask) is shown in 
the lower right with the visibility splodges corresponding to each 
baseline circled. 



The advantages of this approach are severalfold. By 
preserving non-redundancy, a given baseline in the pupil 
formed by any pair of subapertures translates uniquely 
to a single spatial frequency in the detector plane. The 
fringe observables are extracted from the Fourier trans- 
form of the interferogram (see Figure [2J. Each observ- 
able is a fringe complex visibility. The preservation of 
the non-redundancy relation ensures that the extracted 
fringjes can be used to form closure phases ijJennisonl 
Il958t l(^ornwel]llT989l) . The compelling advantage of the 
use of closure phase data is the rejection of any residual 
pupil-plane phase errors, which are the source of both AO 
PSF calibration difficulties and speckle noise. The non- 
rcdundant masking technique therefore rejects the phase 
noise associated with both the instantaneous and time- 



averaged AO system performance. Images can be recon- 
structed using self-calibration techniques. Although the 
dynamic range achieved here is modest by comparison 
with conventional AO it is uniquely close to the central 
star, within a few X/D, which is an area not accessible 
to coronagraphs. The use of closure relations in radio 
interferometry has enabled imaging with d ynamic range 
exceeding 10 5 (e.g. lHardcastle et alJ l|2003j) ). 

2.2. GJ 802 observations 

GJ 802 was observed at the Palomar 200" telescope 
with conventional AO imaging and 9 hole aperture mask- 
ing interferometry on 2004 September 2 UT, in good see- 
ing. Uncompensated images earlier in the night showed 
0'.'6 FWHM seeing at H band (0"75 in V band). The 
conventional AO imaging placed an upper limit of a 
contrast ratio of 0.0 5 for any companion at ^100 mas 
l)Pravdo et al.l I2005D . Imaging observations have also 
been attempted with Keck/LGS (C. Gelino, pers. comm) 
and HST (GO-10517) but have not detected the compan- 
ion. 

Interferograms were recorded using the Fowler sam- 
pling mode of PHARO on a 256x256 subarray. PHARO 
provides a mode whereby all reads of a Fowler sampling 
sequence can be saved. We use this mode and the mini- 
mum exposure time to save a data cube of sixteen sequen- 
tial non-destructive reads of the detector without reset. 
This provides fifteen pairwise 431 ms exposures in a total 
exposure time of 6465 ms, thus very efficiently recording 
a large number of short exposures, so long as the detector 
does not saturate in the total exposure time. For these 
observations, we recorded 115 image cubes yielding 1725 
431 ms exposures on source. The large number of frames 
allows good estimation of the errors. The observations 
show 0.6° RMS closure phase scatter. 

Calibration of the interferograms is achieved by observ- 
ing an unresolved source to measure the system visibili- 
ties. It is usually considered necessary to choose a cali- 
brator that is as similar as possible to the target star in 
both wavefront sensor (approximately R band) and sci- 
ence camera (H band) brightness, at similar airmass. We 
select calibrator stars by sea rching the USNO CCD As- 
trograph Catalog (UCAC2) ijZacharias et alJl2004|) and 
2MASS catalog for stars nearby in the sky with similar 
properties. The UCAC2 catalog bandpass is between V 
and R bands, and for practical purposes has proven to be 
a similar magnitude scale to the PALAO wavefront sen- 
sor. For these observations 2MASS 20494024+4526398 
(2UCAC 47204238; UC=14.72 mag; H=7.73 mag) was 
selected as a calibrator star. 

2.3. Data Analysis 

The data was dark subtracted, flatfielded, and anal- 
ysed with a custom software pipeline written in IDL. 
The pipeline outp uts a bispectrum in OIFITS format 
l|Pauls et al.ll20 05). A binary model is fit to the bispec- 
trum with a reduced x 2 method. In practice we have 
found that the visibility amplitude calibration is poor, 
and superior results are achieved with a fit to the closure 
phase alone. Presumably this is because the visibility 
amplitude calibration is susceptible to the same fluctua- 
tions in seeing and AO performance between source and 
calibrator that plague conventional imaging with AO. 
As discussed in section |2~D the closure phase rejects the 
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Fig. 3. — Likelihood function cross-section for GJ 802 binary 
model fit. The contours represent 90, 99, and 99.9% confidence 
levels. 



TABLE 1 
Orbital Elements 



Quantity 



Value 



Absolute Parallax 
Proper Motion 
Position Angle 
Period 
Total Mass 
Semi-major Axis 
Eccentricity 
Inclination 
Long. Of Asc. Node 
Arg. Of Periastron 
Epoch (y) 
Primary Mass 
Secondary Mass 



64.5 ± 2 mas 
1933 ± 1 mas 
26.6° ±0.5° 
3.13 ± 0.04 y 
0.24 ± 0.05 Mq 
1.32 ± 0.09 AU 
0.60 ± 0.29 
82° ± 2° 
18.5° ±4.5° 
222.5° ±25.5° 
2000.51 ± 0.19 
0.175 ± 0.021 Mq 
0.064 ± 0.032 M 



residual phase errors and is therefore expected to be ro- 
bust. Although in principle the closure phase is self- 
calibrating, there are systematic non-zero closure phase 
errors of a few degrees. Therefore it remains necessary to 
calibrate the non-zero closure phases. The source of these 
non-zero closure phases is not entirely understood, but 
the level is consistent with the expected telescope and 
AO system residual wavefront errors and detector flat- 
fielding errors. Once the visibility amplitude is rejected 
from the analysis, it is also possible to include additional 
calibrators observed throughout the night. These addi- 
tional calibrators usually improve the estimation of er- 
rors, and provide a robustness against possibility that 
the calibrator itself is an unknown binary. We include 
observations of HD 4915 (GOV, V=6.76, H= 5.416) HD 
28005 (GOV, V= 6.71, H=5.506) as additional calibrators 
in this analysis. 

Likelihood contours for the binary model parameters 
are shown in Figure The derived binary model is sep- 
aration 102 ± 7 mas at position angle 36.1 ± 4.5 degrees. 
The H band contrast ratio is 74.4 +/- 18.5 (AH = 4.68 
±0.28 mag). 

3. ORBITAL PARAMETERS 

The measured position of GJ 802B, show n in Figure IH 
is very near the expected position from iPravdo et alJ 
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FlG. 4.— Orbit of GJ 802B. STEPS astrometric observations are 
diamonds. The uncertainties shown are the photocenter astrom- 
etry uncertainties multiplied by the ratio (~4) of the displayed 
Keplerian orbit (curve) to the photometric orbit. The resolved AO 
observation is the triangle. 



(|2005J) although this prior information is not involved 
in the fit to the closure phases used to extract the as- 
trometry. Updated orbital elements combining this re- 
solved obse rvation and addition al astrometric observa- 
tions since IPravdo et all ffiiol are shown in Table ^ 
This orbit is consistent with the pure photo center as- 
trometry orbit derived in IPravdo et alJ l)2005f) . 

4. DISCUSSION 

The luminosity of GJ 802B can be determined pre- 
cisely by a differ ential measuremen t from GJ 802A. The 
2MASS catalog l|Cutri et alJl2003j) records the H band 
brightness of GJ 802 as 9 .058± 0.019 mag. A dopting 
the parallax determined bv IPravdo et al.l l)2005j) of 64±2 
mas, the absolute magnitude of GJ 802 is M# = 8.11 
± 0.07 mag. Using the H band contrast ratio of AH = 
4.68 ± 0.28 mag, we determine the absolute magnitude 
of GJ 802 B to be M H = 12.79 ± 0.3. Comparison with 
models of IBaraffe et alJ l|2003j) admits a large range of 
possible masses depending on age (see Figure |3J), 0.057- 
0.074M Q , with Teff = 1850±50K for models of ages 1-10 
Gyr. With an age estimate of < 6 Gyr based on activity 
(Pravdo et al. 2005), the Baraffe models indicate a mass 

- 0.07 M . 

For models with age > 5 Gyr, the mass range con- 
sistent with this luminosity is remarkably narrow, 0.072 

- 0.074. This modcl-dcpcndcnt mass range is narrower 
than the present dynamical mass determination s based 
on the STEPS orbit alone (jPravdo et all 12005(1 or this 
work. Although this mass is consistent with the orbital 
solution based on the STEPS and AO masking result, 
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the high mass required by the models if GJ 802 is old 
demands an unusually high eccentricity (e >0.8) for the 
astrometric orbital solution. It is therefore tempting to 
conclude that the GJ 802 is young (< 1 Gyr), which 
would admit a lower eccentricity for the orbital solution. 

A sample of field objects selected with a luminosity 
near the substellar limit would be dominated by stellar 
objects since the cooling time of low mass stars dramat- 
ically exceeds that of brown dwarfs. Further, since the 
cooling time increases with mass, the distribution of field 
objects below the substellar limit contains many more 
old massive brown dwarfs than young low mass objects. 
The conclusion that GJ 802 is young would be remark- 
ably puzzling since then the only three resolved binaries 
with dynamical masses below 0.08 Mp (GJ 802B; GL 56 9 
Bab (Za patero Osorio et al.l 120041: IS imon et all 1 2006): 
2MASSW J0746425+2000321AB lIBouv et al.ll200lv rare 
inferred to be young brown dwarfs, despite the fact that 
young low mass objects are less likely to be found than 
old high mass ones. These conclusions are suggestive 
that the models are under-predicting the luminosity of 
substellar objects. 

Ultimately, further observatieoons will provide model- 
constraining dynamical measurements of the masses of 
both the primary and secondary. If GJ 802 is an old 
star, then GJ 802B is remarkably close to the brown 
dwarf/substellar boundary, and more measurements with 
this technique will provide tight constraints on the sub- 
stellar evolutionary models. 



Fig. 5. — Observed mass a nd luminosity of GJ 802B compared to 
theoretical isochrones from (Baraffe et al. 2003). The mass is the 
dynamical mass on the basis of the orbital fit shown in Table IT1 and 
Figure Uncertainty in the mass is large due to the poorly con- 
strained orbital eccentricity. Future observations should constrain 
the eccentricity and therefore the mass precisely. 
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